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The fundamental theory and the semiclassical description of loop quantum cosmology (LQC) 
have been studied in the Friedmann-Robertson- Walker and Bianchi I models. As an extension 
to include both anisotropy and intrinsic curvature, this paper investigates the cosmological model 
of Kantowski-Sachs spacetime with a free massless scalar field at the level of phenomenological 
dynamics with the LQC discreteness corrections. The LQC corrections are implemented in two 
different improved quantization schemes. In both schemes, the big bang and big crunch singularities 
of the classical solution are resolved and replaced by the big bounces when the area or volume scale 
factor approaches the critical values in the Planck regime measured by the reference of the scalar 
field momentum. Symmetries of scaling are also noted and suggest that the fundamental spatial 
scale (area gap) may give rise to a temporal scale. The bouncing scenarios are in an analogous 
fashion of the Bianchi I model, naturally extending the observations obtained previously. 
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I. INTRODUCTION 

It has been long suggested that the singularities in general relativity signal a breakdown of the 
classical theory and should be resolved by the quantum effects of gravity Loop quantum gravity 
(LQG) is one of such candidate theories of quantum gravity and its application to cosmological 
models is known as loop quantum cosmology (LQC) (see p] for a review). The comprehensive 
formulation for LQC has been constructed in detail in the spatially flat and isotropic model with a 
free massless scalar field [21 HUH], showing that the quantum evolution is deterministic across the deep 
Planck regime and the cosmological singularity is replaced by a big bounce for the states which are 
semiclassical at late times. This construction was then extended to k — ±1 Friedmann-Robertson- 
Walker (FRW) models to include intrinsic curvature [HI E] as well as Bianchi I models to include 
anisotropy El El HD] either in the fundamental theory of LQC or at the level of phenomenological 
dynamics; the studies in extended models affirm the resolution of cosmological singularities and the 
occurrence of big bounces. 1 

To further extend this formulation and enlarge its domain of validity, the next step is to investigate 
loop quantum geometry of the black hole and to see whether the black hole singularity is also 
resolved. The simplest step is to consider the interior of a Schwarzschild black hole, in which the 
temporal and radial coordinates flip roles and thus the metric components are homogeneous with 
the Kantowski-Sachs symmetry. Because of homogeneity, the loop quantization of the Schwarzschild 
interior can be formulated in a similar fashion of LQC. This has been developed in [TH [121 H3] and 
its phenomenological dynamics studied in [H] shows that the black hole interior is extended to a 
white hole interior through the bounce, which resolves the singularity. 



'Electronic address: chiou@gravity.psu.edu 

1 However, not all extended models have been developed fully rigorously and in some models the affirmation is only based on the 
semiclassical treatment, the validity of which has been justified only for the spatially flat and isotropic model with a free massless 
scalar field. With this caveat, it would be more precise to call the heuristic treatment of semiclassical approach used in this paper 
" phenomenological dynamics" instead of "effective dynamics" as it has not been shown to be so. (The author thanks Martin Bojowald 
for this comment.) 
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The analysis in [T3] is based on the original quantization strategy (referred to as the "/i -scheme" 
in this paper) used in |12j . which, as a direct transcription of the original LQC construction in [2], 
introduces a fixed parameter to impose fundamental discreteness of quantum geometry. However, 
it has been argued that the /i -scheme quantization in LQC leads to a wrong semiclassical limit 
in some regimes and should be improved by replacing the discreteness parameters with adaptive 
variables which depend on the scale factors [I] . Two improved strategies for loop quantization of the 
Schwarzschild interior were investigated in jl5j at the level of phenomenological dynamics, revealing 
that the black hole singularity is resolved with the correct semiclassical behaviors in the regime away 
from the black hole singularity and the event horizon. 

However, due to the absence of matter content, the result of the Schwarzschild interior is not 
easy to be directly compared with the bouncing scenario of LQC. In particular, the occurrence of 
bounces of triad variables is not indicated by the matter energy density and difficult to be pinpointed; 
moreover, the interesting symmetry of scaling observed in [TU] is not notable. 

In order to bridge LQC in the isotropic and Bianchi I models with the case of Kantowski-Sahcs 
symmetry, instead of focusing on the Schwarzschild interior, we consider a cosmological model of 
Kantowski-Sachs spacetime by introducing a free massless scalar field. The inclusion of the scalar 
field gives a few advantages. First, in the presence of the scalar field, there is no event horizon 
and therefore we are working on a "self-contained" cosmological model with no need to refer to the 
"exterior" (see Remark|2]in Sec. II B). Second, the scalar field serves as emergent time and makes our 
model directly analogous to those which have been carefully studied. Third, the scaling symmetry 
noted in [TU] can be analyzed in a similar manner with the reference of the matter momentum. 
Furthermore, the Kantowski-Sachs spacetime as a cosmological model possesses both anisotropy 
and intrinsic curvature; in a sense this is a hybrid of the Bianchi I and k = +1 FRW models (see 
Remark [I] in Sec. II B) and thus sets a new testing ground for the LQC ramifications. 

Based on the same semiclassical approach of [10J to incorporate the LQC discreteness corrections, 
the phenomenological dynamics of the cosmological model in Kantowski-Sachs spacetime with a free 
massless scalar field is investigated in this paper with two improved quantization strategies (called 
"/i-scheme" and "/i'-scheme"). The investigation shows that both the big bang and big crunch 
singularities of the classical solution are resolved and replaced by big bounces and as a result the 
evolution of phenomenological dynamics follows (semi)-cyclic patterns. 

As a direct analog of the Bianchi I model in [TU], the /i-scheme and /i'-scheme give rise to different 
bouncing scenarios, distinct from each other not only in detail but also qualitatively. In the Jx- 
scheme, the indications of the occurrence of big bounces are "directional densities" Qb and g c and 
the triad variable pb is perfectly periodic while the other triad variable p c only bounces a few times 
and grows to infinity in the far future and past. By contrast, in the /i'-scheme, it is the matter 
energy density p<f> that signals the big bounces and therefore pb and p c get bounced roughly around 
the same moment. Additionally, the /i'-scheme has the problem that p c eventually descends into 
deep Planck regime (while pb grows huge), signaling a breakdown of the semiclassical description. 

As in the Hamiltonian framework for homogeneous models, we have to restrict the spatial inte- 
gration to a finite sized shell X x S 2 to make the Hamiltonian finite. This prescription raises the 
question whether the resulting dynamics is independent of the choice of the finite interval X. It 
can be shown that the phenomenological dynamics in the /i'-scheme is completely independent of 
the choice of X as is the classical dynamics, while the phenomenological dynamics in the /i-scheme 
reacts to the macroscopic scale introduced by the boundary condition of X. This is an important 
difference between these two schemes (although the independence of X is not required with quantum 
corrections). 

In addition to the issues related to the dependence on X, the phenomenological dynamics also 
reveals interesting symmetries of scaling, which are reminiscent of the relational interpretation of 



3 



quantum mechanics. It is also suggested that the fundamental scale (area gap) imposed for the 
spatial geometry may gives rise to a fundamental scale in temporal measurement. This observation 
further supports the speculations made in [TU] . 

This paper follows the steps in [TU] as closely as possible and uses notations in the similar style. 2 
In Sec. [IT], the Ashtekar variables of Kantowski-Sachs spacetime are introduced and the classical 
dynamics with a massless scalar source is solved in Hamiltonian formalism. The phenomenological 
dynamics with LQC discreteness corrections is constructed and solved in Sec. Ill for the £L- and p,'- 
schemes, respectively. The scaling symmetry and issues about relational measurements are discussed 
in Sec. [TV} Finally, the results are summarized and discussed in Sec. |V} As a comparison to the p,- and 
/i'-schemes, the phenomenological dynamics in the /i -scheme is also included in Appendix [Aj The 
details of heuristic arguments and motivations for the p- and /i'-schemes are given in Appendix iBj 



II. CLASSICAL DYNAMICS 



In this section, we first briefly describe the Ashtekar variables in Kantowski-Sachs spacetime [T2] . 
In the Hamiltonian framework, we then solve the classical solution in terms of Ashtekar variables 
for the Kantowski-Sachs cosmology with a free massless scalar field. 

A. Ashtekar variables in Kantowski-Sachs spacetime 

The metric of homogeneous spacetime with the Kantowski-Sachs symmetry group R x SO (3) is 
given by the line element: 

ds 2 = -dr 2 + g xx {r)dx 2 + g nn (T)dfl 2 

= -N(t) 2 dt 2 + g xx {t)dx 2 + g ee (t)d9 2 + g H {t)d4> 2 , (2.1) 

where r is the proper time, N(t) is the lapse function associated with the arbitrary coordinate time 
t via N(t)dt = dr and dQ 2 represents the unit 2-sphere given in polar coordinates as 

dn 2 = d6 2 + sin 2 6 d(f) 2 . (2.2) 

The topology of the homogeneous spatial slices is £ = R x S 2 , which is coordinatized by x G R, 
9 G [0,tt] and G [0,2tt]. 

As in any homogeneous cosmological models, on the homogeneous spacelike slice E, we can choose 
a fiducial triad field of vectors °e° and a fiducial cotriad field of covectors °uj l a that are left-invariant 
by the action of the Killing fields of S. (Note O e a i °io % h = 5%.) The fiducial 3-metric of £ is given by 
the cotriad °u) l a \ 

°to = K°^«r (2-3) 
In the comoving coordinates (x, 9, 0), we can choose °q a b to have 

°q ab dx a dx b = dx 2 + d9 2 + sin 2 9 dtf, (2.4) 
which gives °q := det °q a b = sin 2 9. 



2 Assiduous readers are encouraged to look at |10| to see the close parallel (but note that the term "effective dynamics" has been rephrased 
as "phenomenological dynamics" in this paper). Also see Remarklllin Sec. Ill Blfor some comments. 
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In connection dynamics, the canonical pair consists of the Ashtekar variables: the densitized triads 
E i(x) and connections A a l (x), which satisfy the canonical relation: 



{Aj(x), E b ,(x')} = 8ttG 7 5) S b a 6 3 (x - x'), (2.5) 

where 7 is the Barbero-Immirzi parameter. In the case that connections and triads admit the 
Kantowski-Sachs symmetry E x 5*0(3), A a l and E j after gauge fixing of the Gauss constraint are 
of the form [12J: 

A = A a l Tidx a = cr^dx + bT2(16 — br\ sin 6 d<f) + T3 cos 8 d(f), (2-6) 
E = E iTid a = p c r 3 sin 9 d x + p b r 2 sin 6 d e - p h T X 8$, (2.7) 

where b, c, p b , p c are functions of time only and Ti = —i<Ji/2 are SU(2) generators satisfying 
[Tj,Tj] = €ij k Tk (with (Tj being the Pauli matrices.) 

The symplectic structure on the symmetry-reduced phase space is given by the complete symplectic 
structure [as in (2.5)] integrated over the finite sized shell X x S 2 : 

Q = [ d 3 x dA a \x) A dE a i(x) = (dc A dp c + 2db A dp b ) , (2.8) 

87TG7 J IxS 2 2G7 V J 

where the integration is over 9 G [0, 7r], <f) G [0, 2n] and restricted to x G X := [0, L}; the finite interval 
X is prescribed to circumvent the problem due to homogeneity that the spatial integration over the 
whole spatial slice R x S 2 diverges. [We will see that this prescription does not change the classical 
dynamics but might have effects on the quantum corrections.] The reduced symplectic form leads 
to the canonical relations for the reduced canonical variables: 

{b,p b } = G 1 L- 1 , {d,p c } = 2G 1 L- 1 (2.9) 

and {b, c} = {pb,Pc} — 0. It is convenient to introduce the rescaled variables: 

b := b, c := Lc, p b := Lp b , p c := p c , (2.10) 

which satisfy the canonical relations: 

{b,p b } = G 1 , {c, Pc } = 2G 7 . (2.11) 

The relation between the densitized triad and the 3-metric is given by 

qq ah = ^E a t E b v (2.12) 

which leads to 

9na = 9ee = 94,4, sin 2 9 = p C) g xx = -y—. (2.13) 

Let S X( p, S x g and Sg^ be the three surfaces of interest, respectively, bounded by the interval X and 
the equator, X and a great circle along a longitude, and the equator and a longitude (so that Sg^ 
forms a quarter of the sphere S 2 ). It follows that the physical areas of S x ^, S x $ and Sg^ are given by 

= A x9 = 2nLy/g xx g nn = 2np b , A d(t) = ng nn = np c , (2.14) 

and the physical volume of X x S 2 is 

V = inLy/g^gnn = A-np b ^/p~ c . (2.15) 
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This gives the physical meanings of the triad variables pb and p c . 3 

The gravitational part of the Hamiltonian constraint of Kantowski-Sachs spacetime is given in 
terms of Ashtekar variables as 



H, 



N 



grav 



2G 7 2 



2bcy/p~ c + {b 2 + Y 



Pb 

\fPc 



(2.16) 



This can be derived from the Hamiltonian constraint of the full (unreduced) theory. (See the text 



toward (B9) in Appendix [B] or [T2"] 



B. Classical solution 

The vacuum solution of Kantowski-Sachs spacetime is identified with the interior of the 
Schwarzschild black hole [I2]- 4 In this paper, in order to extend the results studied in pU] for 
the Bianchi I model to the Kantowski-Sachs spacetime, instead of the vacuum solution, we inves- 
tigate the cosmological model of Kantowski-Sachs spacetime with the inclusion of a homogeneous 
massless scalar field <f)(x,t) = <fi(t) without introducing any potential of <fi (i.e. <fi is free). 5 

The total Hamiltonian constraint is given by the gravitational part H gvav plus the scalar field part 
H^: that is 



H = H gTSLV + Hcf, = ifgrav + NYp^ 



N 



2G 7 2 



2bcVp c + (b 2 + 1 2 )^ 



where the matter energy density ps is given by 



2V 2 32n 2 p 2 p c 



with ps being the scalar field momentum: 



Np 2 



(2.17) 



(2.18) 



V0 := V 



dr 



(2.19) 



To solve the classical solution, we can simplify the Hamiltonian by choosing the lapse function 
N = lPb\fp~c corresponding to the new time variable dt' = (7Pby / p^) _1 rfT. The rescaled Hamiltonian 
is given by 



H' 



1 



2G 7 



[2bcp bPc + (b 2 + 7 2 )p fe 2 ] + 7 



In 



(2.20) 



3 More precisely, in l |2.14| l pi, and p c should be |p(,| and \p c \ |12| . With the gauge fixing pi, > 0, the opposite sign of p c corresponds to the 
inverse spatial orientation, which we do not need to consider in this paper. 



In the standard Schwarzschild solution, the metric of the black hole interior is of the form of j2.l| : 



ds 2 



-N{tfdt 2 + 9xx (t)dx 2 + g nQ {t)dft 2 



/2GM 



dt 2 



(2GM 



1 )dx 2 + t 2 dQ 2 



where t (£ [0, 2GM], x £ R and M is the mass of the black hole. The black hole singularity corresponds to t = and the event horizon 
corresponds to t = 2GM. 



Do not confuse the matter field <j>{i) with the polar coordinate <f> in (2.1 1 



The equations of motion are governed by the Hamilton's equations: 

> Pcf, is constant 

{<P,H} 



dp,!, 
df 



{P<P,H} = 

7 



df 
dc 

df 
dp c 
df 
db 

~df 
dpb 
df 



Air 



P<t>, 



{c,H'} = 2G 1 



dH' 

dp c 



-2cbpb 



d H' 

-2G*7 -s- = ZPcbpb, 



dc 



{b,H'} = G 1 



dH' 



-b (bp b + cp c ) - j 2 p b , 



dp b 
d H 1 

{p b , H'} = -G7 -jy- = p b (bp b + cp c 



as well as the constraint that the Hamiltonian must vanish: 

Gi 2 pI 



H' = 



47T 



2bcp b pc + {b 2 + 1 2 ) p 2 b . 



Notice that substituting ( |2.14[ ) into ( |2.24[ ) and (]2.26[) gives us 
b 



7 dp c d 

o i—~r~ = J~rv9nn 
2y/p c CLT CLT 

7 dp b ^Pb dp c 



V2 dr 



P 



, 3 / 2 dr 



2p 



7 



d f Pb 



dr 



7^ (^Vfe) 



(2.21) 
(2.22) 
(2.23) 
(2.24) 
(2.25) 
(2.26) 

(2.27) 

(2.28) 
(2.29) 



which tells that, classically, the connection variable b is the time change rate of square root of the 
physical area of S 2 (up to constant (47r) -1 7) and c is the time change rate of the physical length of 
I (up to constant 7). 



To solve the equations of motion, first note that combining (2.23) and (2.24) gives 

d 



dt 



y(Pcc) = p c c = ^K c is constant, 



and on the other hand, (2.25) and (2.26) yield 

d 



dt 



7 (K b 



-jpl with p b b =: jK b (f). 



The Hamiltonian constraint (2.27) then reads as 

Kl = 2K b K c + K 2 +p 2 b 

if we define 

Gpl 



By (2.31) and (2.32), we have 



7" 



-i_ dK b 
df 



Kj 



Ka dKh 



Air 



2K b K c + K 2 - Kl 



(2.30) 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



where (2.22) has been used. It follows from (2.22) that is a monotonic function of time and 



therefore can be regarded as emergent time. In terms of 0, the solution to (2.34) is given by 



K b {<j>) = -K c - JKl + K\ tanh 



4ttG {Kl + K}) 
Kl 



-(po)-a 



with a being the constant specified by the initial state: 



a := tanh" 



_! / K h {(j> Q ) + K c 



Kl + Kl 



In terms of K b ((p) and the constants K c and K^, (2.24) and (2.26) read as 



1 dp c V^G 
— -TT = 2 K b (<p), 
p c d(p K<p 



1 dp b _ v / 4ttG 
p b Ka, 



[KM) + # c 



the solutions to which are given by 
Pb{4>) = Pb(4>o) cosh a 



cosh 



4ttG [Kl + XJ) 
Kl 



-<po)-a 



and 



P c (0) = 9an((f>) 



4VnGK c 



Pc(4>o) (cosh a) e K ^ 



cosh 



Consequently, we have 



9xx{4>) 



Pb Pb((f>of 



4nG {Kl + Kf) 
Kl 



■<<t>-M 



i ) - « 



L 2 p c L 2 p c ((J) ) 



and 



V(0) = 47rp 6 (0 o )v / P c (0o) (cosh a) 2 e ^ 



2VnGK c 



2 _~ 'itT'i Hl)) 



X 



cosh 



4vrG (Kl + Jg) 



Jo) - a 



Furthermore, from (2.35) and (2.39), we have 

-Kr~ 



K b {<j>) 



K b1 



Kl + Kl-pl for0> 



K c +JKl + Kl-pl for0< 



(2.35) 



(2.36) 



(2.37) 
(2.38) 



(2.39) 



(2.40) 



-2 



(2.41) 



(2.42) 



(2.43) 
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with the convention K$ > 0. 6 That is, Kf, 
Kb = in the expanding phase (i.e. < m 



K^ in the collapsing phase (i.e. 
), while Pb reaches the maximal value 



> 



Pb, 



and 



(2.44) 



at the epoch = ma; 
equation of K h given by (|2.32|). 
K b ±: 



Note that the two branches K^i are the two solutions of the quadratic 
In particular, when p\ <C K 2 + K?, K b approaches the constants 



K h 



K b+ {K C 
K b -(K C 



-K e -Jia + K* 



-K e +jKZ + K* 



as > max , 
as < max , 



(2.45) 



and we call the periods with pi -C K\ + K\ the "if^t-asymptotic phases". 7 
Additionally, the stress tensor of the scalar field is given by T ab = V a 0Vfe0 • 
the Einstein's equation R a b — \Rg a b = SnGT a b gives the scalar curvature 



R = -8nGT = -8nGT a a = 8ttG0 2 = 8nG 



V 2 ' 



^abV c 0V c and thus 



(2.46) 



In the future and past, the classical solution eventually approaches singularities with the asymp- 
totic behaviors: 



Pb -> 0, p c 



and 



gnn 



9.i- 



0, V^O, R 



oo 




as K c 
as K r 



oo, 
— oo. 



as 



±oo 



(2.47) 



(2.48) 



The finite sized shell X x S 2 collapses to a point at K c cf) — > — oo while to an infinite line at if c — > oo. 
Regardless the different signatures, we always call the singularity in the future (0 — > oo) the 
crunch singularity and the other in the past (0 — > — oo) the bang singularity. Although the 
singularities correspond to ±oo in 0, the universe actually takes finite proper time to reach either 



of the singularities (in forward or backward evolution), which can be verified by showing that (2.49) 
is finite even when the upper limit of the integral is taken to be = ±oo. The behaviors of the 
classical solution are depicted in Fig. [TJ 

It should be noted that, if we flip the sign of K c (while fix K^), we have Kb(4>) — ► —Kb(—<f)). 
Equations (2.39) and (2.40) then tell us that K c — ► — K c corresponds to the time reversal. 

Notice that, by (2.14), (2.19) (2.28) and (2.29), p^, p b and c depend on the choice of the interval 
X and scale as p^,pb,c oc L, while p c and b are independent of X. As a result, the constants of 
motion K^ and K c (as well as the function Kb) all scale as oc L. The ratios Kc/K^ and Kb/K^ are 
nevertheless independent of X and hence the classical solutions of Pb/Pb{(f>o) and p c /Pc(4>o) given by 
(2.39) and (2.40) do not depend on the choice of X. Furthermore, once p b (<p) and Pb{4>) are solved, 
to know the solutions Pb(r) and P&(t) as functions of r, we only need to convert back to r via 



T 



T~0 



dr 



<p y 
o V0 



^PbJp~ c 



(2.49) 



6 We get the equivalent solution if the signs of K^, K c and are flipped simultaneously. Thus, we fix the convention > 0. 

7 We will show that the loop quantum corrections (in both the p.- and ft! '-schemes) take effect only in the iff,±-asymptotic phases, provided 
that the semiclassicality condition is well retained. See footnotes 1121 If 81 and If 91 
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(a) (b) 




-10 12 -2-1012 



FIG. 1: Classical solution. The initial condition is given at (f>o = <^max with pb(4>o) = Pt>,max = y K£ + K$ = 5. x 10 5 4>i 

(£pi := VGh), p c (0o) = 4. x 10 5 4i, K c = 3. x 10 5 4i and JQ, = 4. x 10 5 4i (i.e. p</> = 4. x 10 5 /iV4^G). (a) p b (0) and 
p c (<f>) = gnsi{<t>)- (b) g xx ((j>)- (c) V(0). (d) Kb(4>), which asymptotically approaches iff,± as indicated by dashed lines. 



where, again, the dependence of L is gone. Therefore, whether in terms of the proper time r or the 
emergent time 4>, the classical dynamics is completely independent of the finite interval X we choose 
to make sense of the Hamiltonian formalism. [The independence of the choice of X is not necessarily 
retained when quantum corrections are taken into account. In any case, however, the dynamics is 
independent of X in terms of r if and only if it is so in terms of 0.] 

Remark 1 Comparison with the Bianchi I model 

It is instructive to note that the Hamiltonian formulation with the Kantowski-Sachs symmetry 
features some characteristics reminiscent of the Bianchi I model. Exploiting these resemblances 
can help us solve and understand the equations of motions for both classical and phenomenological 
dynamics. 

Recall that in [TUI the classical Hamiltonian constraint in the Bianchi I model is given by 



H 



C2P2C3P3 + C1P1C3P3 + C1P1C2P2 



+ 



(2.50) 



If we formally make the identifications: 

Pl,P2 - 
P3 - 



Pc, 



Cl,C 2 

c 3 — 



b. 



Aire, 



(2.51; 



then (2.50) becomes the same as (2.11) (with N = 1), except that the term ^PbPc 1 ^ 2 inside the 



bracket of (2.17) is missing. Since the term ^PbPc 1 ^ 2 corresponds to the intrinsic curvature of S 2 
(see Appendix \B\), it suggests that if we ignore the curvature of S 2 , the Kantowski-Sachs spacetime 
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is the same as the Bianchi I model as I = 1,2 directions are identified with the spherical directions 
8q and and 1 = 3 with the "radial" direction d x . 

The presence of the curvature of S 2 complicates the equations of motion and makes the evolution 
look somewhat familiar with the k = +1 FEW model, which has both expanding and collapsing phases 
and also possesses both the big bang and big crunch. Nevertheless, in the asymptotic phases, 
the S 2 curvature is negligible and the strategy used to solve the dynamics of Bianchi I models as in 
pm can be carried over. In particular, we exploit the fact that in classical dynamics, p c c is constant 
and so is pbb approximately in the Kb±- asymptotic phases. 8 Since the quantum effect is expected 
to take effect only in the Kb±- asymptotic phases (in which pb is small enough), this strategy can 
be easily adopted to deal with the phenomenological dynamics in both the fx- and JjC ' -schemes in the 



same fashion as in [10]. In Sec. Ill, we will follow closely the treatment used in UDf to analyze the 
phenomenological dynamics. 

It should also be noted that, as discussed in ^ for JCi in the Bianchi I model, K c and Kb char- 
acterize the anisotropy (between the spherical directions and the radial direction.) The Hamiltonian 



constraint (2.32) can be understood as the relation which relates matter energy with anisotropy and 



spherical curvature. 

Remark 2 Comparison with the Schwarzschild interior 

If we did not include any matter content, the vacuum solution with Kantowski-Sachs symmetry 
would represent the interior of a Schwarzschild black hole JT^, which encounters the black hole 
singularity (with Pb,p c 0) in the future and the event horizon (with pb — > 0, p c — > (2GM) 2 ) in the 
past. The horizon is however not a singularity but is to be extended to the exterior of the black hole; 
thus, there might be no quantum corrections near the event horizon. (The analysis in J75) / shows that 
the quantum corrections for the event horizon are present only in the p! ' -scheme phenomenological 
dynamics but absent in the fl-scheme.) 

With the presence of a free massless scalar field, the situation is a bit different. The Kantowski- 
Sachs spacetime with a massless scalar possesses the big bang singularity in the past and the big 
crunch singularity in the future but no event horizon anywhere. Therefore, this is a "self-contained" 
cosmological model and the appeal to the "exterior" is not needed. 9 Both singularities are genuine 
singularities (pb,Pc — > 0) an d consequently we shall expect both of them are resolved and replaced 
by the big bounces if the loop quantum corrections are taken into account. However, these two 
singularities have different signatures: one of them yields g xx — > oo while the other gives g xx — > 0. 



III. PHENOMENOLOGICAL DYNAMICS WITH LQC DISCRETENESS CORRECTIONS 

In the fundamental loop quantum theory of Kantowski-Sachs spacetime, the connection variables 
b and c do not exist and should be replaced by holonomies [12J. At the level of phenomenological 
theory, to reflect the quantum corrections on the states which are semiclassical when the universe is 
large, following the procedures used for the isotropic cosmology [16] and the Bianchi I model [10], 
we take the prescription to replace b, c with 

sin(/2 6 b) sin(/x c c) 
b — ► — z , c — ► — z , (3.1) 

H He 



8 Cf. The slightly different notations are used in 1101 . in which cjpj =: SvrG'yhK.j and p$ =: hy/SnG/Cgy. The constants K.j = ICkj and 
K,$ = K.n$ are to be identified with Kb, K c and K,p used in this paper up to overall dimensionful constants. 

9 The reader might think the Kantowski-Sachs spacetime with inclusion of matter describes the interior of a collapsing black hole. This 
is not the case, since what we include is a homogeneous matter filed, whereas the collapsing black hold requires inhomogeneous (but 
spherically symmetric) distribution of matter density. 
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introducing the variables fib and fi c to impose the fundamental discreteness of quantum geometry. 10 
The heuristic argument starting from the full theory of LQG for this prescription is presented in 
Appendix [B] 

The discreteness of quantum geometry also modifies the cotriad component oj c := Pbl \fp~c = 
Ly/g xx . The eigenvalues of Cj c are very close to the classical expectations far away from the classical 
singularities but become significantly different from the classical values close to the singularity at 
which pb/ ^Jp~ c diverges [12]. In the semiclassical description, however, this modification on the cotriad 
io c is less important and we will ignore it by simply taking the classical function pb/ ^Jp~ c for the cotriad 
u c . [We will see that for the solutions which are semiclassical far away from the singularities, the big 
bounces take place when p c is still much larger than the square of Planck length and the discreteness 
correction on u c is yet to be considerable. It is the "nonlocality" effect (i.e., using the holonomies) 
that accounts for the occurrence of big bounces.] 11 



As a result, with the prescription of (3.1) adopted and the cotriad component u o c unc hanged, by 
choosing N = r )Pb\fp~ c and dt' = (^Pby/Pc)~ l dT, the (rescaled) classical Hamiltonian (2.20) is modified 
to serve as the effective Hamiltonian for the semiclassical theory: 

= 1 | 2 sin( fe 6) sinfeo) / sin^V | pj 

* 2G7 [ fib fic V fib ) j 8tt 

The phenomenological theory prescribed here is only heuristic and questionable; a more rigorous 
understanding of the quantum dynamics would require more sophisticated refinement. Nevertheless, 
the fact that this phenomenological theory could provide an accurate approximation (for the case 
that the back-reaction is negligible) has been evidenced in the isotropic cosmology (H [161 H3 HH] 
and also affirmed in the Bianchi I model |19j . 

As for imposing the fundamental discreteness of LQG on the formulation of LQC, the original 
construction (/z -scheme) is to take fib and fi c as constants (referred to as fi%,fi° in Appendix |A| and 
5 in [121 HH HSJ)- However, it has been shown in both isotropic and Bianchi I models that the 
/i -scheme can lead to the wrong semiclassical limit and should be improved by a more sophisticated 
construction (/i-scheme) in which the value of discreteness parameters depends adaptively on the 
scale factors (e.g. fi oc l/^/p is used in [1]) and thus implements the underlying physics of quantum 
geometry of LQG more directly [U [10] . 

To impose the discreteness in the Kantowski-Sachs spacetime, there is a variety of possibilities 
in the improved (fi-) scheme. Among them, two well-motivated constructions (referred to as the 
"/i-scheme" and "/i'-scheme" ) are focused in this paper: 

• /i-scheme: 

[a _ [a 

fib = \ — , fic = \ — , (3-3) 

V Pb y Pc 

• fi'- scheme: 

' fa AS 

fib=\ -i ^ = ^~- (3 - 4) 

V Pc Pb 



This prescription is sometimes referre d to as "polymerization" or "holonomization" in the literature. 

1 However, as will be seen in Sec. |III 5] in the /i'-scheme, p c eventually descends into the deep Planck regime at some point and there we 
can no longer trust the phenomenological theory without taking into account the modification on the cotriad u) c . 
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Here A = 2 (^^pi) * s ^ e area 9 a P * n ^ e theory of LQG with £p\ = V Gh being the Planck 
length. 

Either scheme of them has its own advantages and disadvantages and until more detailed physics 
is investigated it is arguable which one makes better sense. In particular, the /2-scheme (in the ver- 
sion for the Bianchi I model) is suggested in [7], since in the construction of the fundamental LQC 
the Hamiltonian constraint in the /2-scheme gives a difference equation in terms of affine variables 
and therefore the well-developed framework of the isotropic LQC can be straightforwardly adopted. 
(However, it is argued in [20] that the /2-scheme may lead to an unstable difference equation.) By 
contrast, the ^'-scheme does not admit the required affine variables and the fundamental LQC of it 
is difficult to construct. On the other hand, the /i'-scheme has the virtue that its phenomenological 
dynamics is independent of the choice of X as will be seen (although this virtue is not necessarily 
required when quantum corrections are taken into account). (More details of the heuristic argu- 
ments and motivations for both schemes as well as their comparison are presented in Appendix [Bj) 
To explore their virtues and differences, we study both the /2-scheme and /Z'-scheme in the con- 
text of phenomenological dynamics in Sec. Ill A and Sec. Ill B respectively. (For comparison, the 
phenomenological dynamics in the /i D -scheme is also presented in Appendix |Aj where we see the 
insensible behavior that pt,, p c at the bounce can be made arbitrarily big.) 



A. Phenomenological dynamics in the (j-scheme 



The phenomenological dynamics in the /2-scheme is specified by the Hamiltonian (3.2) with /2b, /2 C 
given by (3.3). Again, the equations of motion are governed by the Hamilton's equations and the 
constraint that the Hamiltonian must vanish; these are 

{p^, H'^} = =>- p^ is constant, 
{</>> H 'fi} = ^> 

am 

-2 



dt' 



dt' 
dc 

dl' 

dPc 
dt' 
db 
dt' 

dpb 
dt' 



3sin(/i c c) ccos(/i c c) 



{p c ,H'} = -2G 1 



OH' 



dc 



2/i c 



2p c cos (/t c c) 



2 

sin(/2{,6) 
fib 



sin(/i b &) 
fib 



~Pb 



~Pb 



f) H' 

{6,^} = G 7 -^ 

OPb 



3sin(/ife6) bcos(fibb) 



sin(/i b fr) sin(/i c c) 

- Pb ~r Z Pc 
/ife /ic 



-1 2 Pb, 



{Pb,H'} 



dH' 

-G7 -7^- = p b cos(/z fe 6) 



sin(/i b 6) sin(/x c c) 
Pb H = Pc 



as well as 



H', = 



An 



sin(/ifefe) sin(/i c c) 



Hb 



/'c 



PbPc + 



sm(fi b b) 
fib 



Pi 



(3.5) 
(3.6) 
(3.7) 

(3.8) 

(3.9) 
(3.10) 

(3.11) 



[Note that in the classical limit /!&&, /t c c — > 0, we have sm(fibb)/fib — ► b, sin(/t c c)//i c — > c and 
cos(/2;,6), cos(/2 c c) — > 1. By inspection, it follows that (3.7)— (3. 11 ) reduce to their classical counter- 
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parts (2.23)-(2.27) in the classical limit.] Also notice that (3.8) and (3.10) lead to 

sm(ft b b) 1 7 dp c 



sin(/2 e c) 

fie 



cos(/2 c c) 2 P l /2 dr ' 

1 7 dPb 

cos(fi b b) V X J 2 dr cos(/2 c c) 2p 3 J 2 dr 



1 



IPb dp c 



which are the modifications of (|2.28|) and (|2.29|) with quantum corrections. 

sin(/2 c c) 



Combining (3.7) and (3.8), we have 



3 sin(/i c c) 
2/2 c 



ccos(/x c c) 
2 



7 +p c cos(/i c c' 



eft' 



d 
dt' 



Pc 



/Jl c 



0. 



which, in accordance with the classical counterpart (2.30), yields the constant of motion: 

sm(p c c) 



Pe- 



ll,, 



lK c . 



Similarly, (3.9) and (3.10) lead to 



3sin(/i fe 6) 6cos( / u b 6) 
2jh 2 



dp b db 
7 +p b cos{n b b) 



dt 



dt' 



d 

~dt' 



Pb~ 



sin(/i 6 6) 



fib 



In accordance with the classical counterpart (2.31), we define 

sin(/2 6 6) - 

Pfe z =: lK b [t ). 

fib 



The Hamiltonian constraint (3.11) now read as 

Kl = 2K b K c + K 2 b +pl 



and K b satisfies the differential equation: 



7 



. x dK h 
dt 1 



Km dK h 



cos(fi b b) (2K b K c + Kl - K%) . 



Substituting (3.15) and (3.17) into (3.8) and (3.10) yields 



1 dp c 
Pc dej) 

1 dp b 
Pb dej) 



2—- — cos(/i c c)K 6 , 
J\<t> 



cos(fi b b) [K b + K c ] . 



(3.12) 
(3.13) 



(3.14) 



(3.15) 



-7 p b cos(/i 6 6). (3.16) 



(3.17) 



(3.18) 



(3.19) 



(3.20) 
(3.21) 



Here, as in the classical dynamics, <fi is regarded as the emergent time via (3.6). [Also note that, 
as in the classical dynamics, it follows from (3.19) that the flipping K c — > — K c g ives r ise to 
K b {(f>) — > —K b (—(f>) and thus corresponds to the time reversal according to (3.20) and (3.21).] 

Equations (3.20) and (3.21) are the modifications of their classical counterparts (2.37) and (2.38). 
Notice that the presence of the cos(- • • ) terms gives rise to the repulsive behavior of gravity as the 
classical solution approaches singularities at Planckian energy density. More precisely, in the Jl- 
scheme phenomenological dynamics, p c and p b get bounced whenever cos(/2 c c) or cos(/2&6) flips signs, 
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respectively. To find out the exact moment of occurrence of the bounces, we investigate cos(p c c) 
and cos(fibb) in more detail. 
First, by ( |3.15 ), we have 



cos 



(p c c) = ± [l - sin 2 p c c] 1/2 = ± 



1 - 



Pc 



1/2 



± 



1 - 



Qc 



-,1/2 



Qc, crit 



where we define the directional density for the "p c -direction" as 

and its critical value is given by the Planckian density ppi times a numerical factor K^/K 2 : 

Qc, crit •— PPb 



(3.22) 



(3.23) 



with 



Ppi := (87rG?7 2 A) _1 . 



(3.24) 



(3.25) 



[Note that the directional density is of the same dimension as the matter density p^ := p\j (327r 2 p 2 p c 
and thus the name.] Therefore, the bounce in p c occurs whenever q c approaches f? C)C rit- 
Similarly, (3.17) gives 



> 1 /2 

cos(pbb) — ± [l — sin 2 p b b\ = ± 



1 - 



Pi . 



(3.26) 



Assuming that the bounce in p b takes place in the if&±-asymptotic phases (i.e. p 2 <C K 2 + K 2 ) so 
that K b = -K c ± [K 2 C +K}- pi] 1/2 w K 6t := -tf c ± [if 2 + K 2 ] 1/2 (see footnote 

T 1/2 



1 



7 2 ^ 2 ± A^/ 2 



pI 



± 



i 



cos(pbb) ps ± 

where the directional density for the "pft-direction" is defined as 

06 



£>fct,crit 



12 ), we then have 
(3.27) 



4 



327T 2 p 3 b 



(3.28) 



and its critical values are given by the Planckian density pp\ times numerical factors K 2 ± /K 2 : 



£?b±,crit : : 



Ppi- 



(3.29) 



Therefore, the bounce in p b occurs whenever g b approaches £>ft+, C rit in the forward evolution (i.e. the 
bounce resolves the big crunch singularity) and £>&-, C rit i n the backward evolution (i.e. the bounce 
resolves the big bang singularity). 12 



The big bounce in occurs as ~ ^K^A, which in turn justifies the approximation <C K% + Ki we have used, provided that 
we have the scmiclassical condition ~ K c ~ K^-j- 2> 7 2 A. 
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In summary, both the big bang and big crunch singularities are replaced by the big bounces, which 
take place in both p c and pb, whenever g c or g b approaches its critical values at Planckian energy 
density (and thus cos(/2 c c) or cos(/4&) flips signs in (3.20) and (3.21) respectively). Furthermore, 
the differential equations (3.19) and (3.21) are independent of p c and c (the dependence on p c and 



c is only through the constant K c ). [Also note that, with the cos(/^6) term in (3.19), Kb becomes 



flat (dKb/d(j) = 0) exactly at the same time when p b get bounced.] As a result, the evolution of pb is 
unaffected by the varying of p c and is expected to be periodic (with respect to 0): i.e. each "cycle" 
of classical evolution of p b is connected through the quantum bridge with the next /previous classical 
cycle. 13 [On the other hand, the numerical analysis shows that the big bounce of p c occurs only a 
few times (only 3 times in Fig. [2|.] 

Notice that the constants and K c remain the same throughout the evolution. However, this 
does not mean that the parameters used to parametrize the classical evolutions in different classical 
cycles remain unchanged, since the physical meanings of h and c are altered before and after the 



big bounce according to (3.12) and (3.13). In order to characterize the classical behaviors of the 



universe in different classical periods, we define the "effective K c " as 



effective K r 



7 1 Pc 



7 dp b jp b dp c 



V2 dr 



7 1 cos(fjLbb)p, 



2p 
sin(/2 c c) 

fie 



3 / 2 dr 



+ 7 [cos(/i 6 6) - cos(/i c c)] p b 



sin(/i 6 b) 
fib 



= cos(fibb)K c + [cos(/Xft6) — cos(/x c c)] K\ 
and similarly the "effective Kb as 

7 dp, 



effective K b := 7 Pb 



2p, 



V2 dr 



-1 t — \ sin WO /_ sf> 
7 cos{n c c)pb — z = cos(ii c c)K b . 



fib 



(3.30) 



(3.31; 



In the classical period, cos(fibb) ~ cos(/i c c) ~ ±1 and consequently the effective K c becomes 

effective K c — > ±K C . (3.32) 

That is, the effective K c in the classical regime can be either K c or — K c . Likewise, the effective Kb 
in the classical period reaches the constant 



effective Ki, 



±K b± (±K c 



(3.33) 



in the -fr^-asymptotic phases. Furthermore, the evolution also admits a new phase for some periods 
of time when cos(/2&&) ~ — cos(/2 c c) ~ ±1. In these "meta-classical" periods, we have 



effective K c — > ± (K c + 2K b ) 



(3.34) 



which is nonconstant. The new phase is "meta-classical" in the sense that p b evolves classically 
and p c is far away from the Planck regime, yet p c does not follow the classical trajectory. [The 
occurrence of the meta-classical phases seems to suggest that the /2-scheme dynamics is problematic, 
giving wrong semiclassical behavior in some periods, but this may not be the case after all since the 
Kantowski-Sachs symmetry could be altered by the bounce if the bounce is of p b alone. We defer 
this issue for further investigation.] 



This periodic behavior is very similar to that of loop quantum cosmology in k = +1 FRW model [5]. Also see Remark^ in Sec. II B 
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For given initial conditions, the equations of motion can be solved numerically. 14 The numerical 
solution is depicted in Fig. [2] Note that the bounces occur at the moments exactly when g b or 
g c approaches their critical values. Also notice that pb and Kb are perfectly periodic while p c only 
bounces a few times and grows up to infinity in the distant future and past. 

It is noteworthy that the directional densities g c and Qb are the indications of the bounces but 
their quantities are not independent of the choice of X, as we know 



Qb 



p\ 


Pc 


ft _ A W 


327T 2 p 3 b 


" 2p b 1 


2A x6 


v\ 


Pb 


A 2 

12 _ A x9 


32vr 2 p3 


" 2p 2 


P 2A 2 



oc L~ 



\ 2 oc L 2 . 



(3.35) 



Therefore, contrary to the classical dynamics, the phenomenological dynamics in the /2-scheme is 
dependent on the choice of the finite sized interval X. Another subtler dependence on X comes 
from the fact that the constant of motion scale as oc L but K c and Kb scale as oc L only 
approximately [see (3.12), ( |3.13[ ), ( |3.15[ ) and (3.17) and notice that ( |3.12[ ) and (3.13) involve the 
quantum modification with cos(- ■ ■ ) terms]. As a result, the ratios Ky K 2 and K^/K b± are slightly 
dependent on X. (Nevertheless, in the classical regime, the dependence on X through these ratios is 
negligible.) 

The problem not to be invariant under different choice of X is absent in the /i'-scheme as will 



be seen in Sec. Ill B However, we should not dismiss the /i-scheme immediately as it is a common 
phenomenon that a quantum system reacts to macroscopic scales introduced by boundary conditions 
(for instance, the well-known "conformal anomaly" as a "soft" breaking of conformal symmetry). If 
we have good physical input to tell what exactly the space is to be enclosed as X (such as in the 
compactified Kantowski-Sachs spacetime in which X = S 1 x S 2 , instead of 1R x S 2 , or in the lattice 
refining model of [20]), the dependence on X could be rather meritorious than problematic and the 
bounce occurrence conditions (g c = £> c ,crit, Qb ~ Qb±, crit) can be understood as: The physical areas 
Ag^ and A x8 = get bounced when any of them undergo the Planck regime (times a numerical 
factor Kl/K"l or K 2 ± /KT) measured by the reference of the scalar field momentum p^. 



B. Phenomenological dynamics in the //-scheme 



The phenomenological dynamics in the /i'-scheme is specified by the Hamiltonian (3.2) with jib, 
fl c replaced by fi' b , pf c given in (3.4). To simplify the equations of motion, we choose a different lapse 
function N = (pby /Pc) ' 1 associated with the new time variable dt" = Pby/p^dr. With the new lapse, 
the Hamiltonian (3.2) is further rescaled to the simpler form: 

1 



HZ 



2G7 2 A 



2 sin(n' b b) sin(/i' c c) + (sin(/x 6 6)) + A 



+ 



(3.36) 



Because |sin(/i' 7 c/)| < 1, the vanishing of the Hamiltonian constraint H'L = immediately implies 



< 



+ 



32n 2 p 2 p c ~ 8vrG7 2 A 8nGp 



< 4p 



pi 



if p c > 7 2 A. 



(3.37) 



The common numerical methods (e.g. Runge-Kutta method) encounter numerical instability at some point if we directly solve the 
coupled differentia l equ ations J3 J1 6J — ( |3 . ]~0| ) . To bypass this problem, which is only a numerical artifact, we solve the reduced coupled 
equations: j3.19| l, ( |3 - 20[ > and (j^^y f° r three variables: K^, p c and pt- The variables b and c can be obtained afterwards via | |3.15| and 
13071. 
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FIG. 2: Solution in the ^-scheme phenomenological dynamics. With the same initial condition as given in Fig. [l] (and 
the Barbero-Immirzi parameter set to 7 = 1). (a) pb(<j>) an d p c (0) = gnn(4>)- Pb is periodic while p c only bounces 3 times and 
grows up (with wiggled motion) toward infinity in the far future and past, (b) g X x{4>)- ( c ) Qb(<f>), which signals the occurrence 
of big bounces of pb when approaching Qb±,ait given by (3.291 and indicated by dashed lines, (d) q c {4>), which signals the 
occurrence of big bounces of p c when approaching f> c , crft given by (3.241 and indicated by the dashed line, (e) cos(flbb): the 
periodic curve; cos(/7 c c): the curve which flips signs only 3 times (for 3 bounces) and eventually remains —1 in the future and 
+1 in the past, (f) Kb(<f>), which becomes flat whenever pb undergoes the bounce, (g) Effective K c with the constants ±A' C 
indicated by dashed lines. (See ( |3.32[ ).) (h) Effective Kb with the constants K b ±(±K c ) indicated by dashed lines. (See ( |3.33| ).) 
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This suggests that the matter density is bounded above and thus the big bounces are expected 
to occur when the matter density approaches Planckian density (provided that p c remains large 
enough) . 

To know the detailed dynamics for each individual pb and p c , in addition to the Hamiltonian 
constraint, we study the Hamilton's equations: 



~d¥ 
d(p 

dJ' 

dc 

dl" 



dp c 
dt" 
db 
dt" 

dp b 
dt" 



Ptf, is constant, 



^p\pc 

{C ,„» } = 2G7 ^E 
cpf e cos(/Z' c c) sin(/46) bpf b cos{p,' b b) [sm{pf b b) + sin(/2 / c c)] j_ _ &m 

a n 



(3.38) 
(3.39) 

(3.40) 



{p ei flJ} = -2G? 7 
{b,H"} = G 1 



^Ap c ^Ap c 

d Hp 2pf c cos(/2' c c) sin(/i / 6 6) 



Airptp 



2„2 ' 



dc 7A 
dH'p cfi' c cos(jj' c c)sm(ji' b b) G^p\ 



dp b 7 A p b ^plPc ' 

dH P> _ /2' 6 cos(/2' 6 6) [sin(/x' 6 6) + sm{fi' c c)} 



{Pb,H'p} = -G7 



56 



7A 



Note that ( |3.41[ ) and (|3_43j) give us 



7 dp c 



cos(^c) 2p y 2 dr 
sin(p! c c) 1 7 dpb 



7Pft dp c 



P'c 



cos{p b b) p y 2 dr cos(^c) 2^ /2 



(3.41) 
(3.42) 
(3.43) 

(3.44) 
(3.45) 



which are the modifications of (2.28) and (2.29) with quantum corrections. 
Inspecting ( 3.40 )-( 3.43), we have 



dt 



% (PcC-Pbb) = — • 

Pc 



(3.46) 



In accordance with the constant K c and the function K b (t') used for classical solutions in (2.30) and 
(2.31), introducing the time-varying function f(t"), we set 

Pc c = 1 (K c + f(t")) 



and 



where K' b satisfies 



J„ dK 'b 



PbPc 



p b b =: 7 {K' b {t") + f{t")) , 
^ dK' b 



dt" 



7 



-pI, 



(3.47) 
(3.48) 

(3.49) 



which is to be compared with the classical counterpart (2.31). To start with a classical regime, we 
set K c « 7 _1 PcC and / ~ 0. 
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Taking (3.47) and (3.48) into (3.36), we have the complicated expression for the Hamiltonian 
0: 



constraint HL 

H 1 



AfKj 

PbPc 



2 sin 



+ sin 




PbPc 



sin 



'A 7 2 

pIpc 



(K c + f) 



pIp 



Pc 



which reduces to 



K 2 = 2K' b K c + K' b 2 + p 2 b 



(3.50) 



(3-51) 



in the classical limit as p\p c 3> A^y 2 K 2 ~ A^y 2 K 2 ~ A^ 2 K b± and with / ~ 0. 

As in the /i-scheme, it is expected that both big bang and big crunch singularities are resolved 
and replaced by big bounces, which bridge one cycle of classical evolution with the next /previous 
classical cycle. One might think that in any classical cycles p c c becomes constant and so should /. 
This is however not necessarily true (but in fact, / ~ only in the particular classical cycle in which 
the initial condition is specified). 15 The reason is that, due to the quantum modifications in (3.44) 



and (3.45), pbb and p c c have different physical meanings before and after the big bounce and thus p c c 



may no longer be constant in the consecutive classical cycle. Instead of p c c, what becomes constant 



in any classical regimes is the "effective K c n . Similar to (3.30) and (3.31), we define 



effective K c := 7 p c 



7 dp b 'jpb dp c 



i 1 ' 2 dr 



7 1 cos(n' b b)p, 



p c ' 2p c 



sin(/x' c c) 



3 / 2 dr 



f'c 



+ 7 [cos(/4&) - cos (ji' c c)] Pb 



siri(// 6 &) 



and 



effective Kb 



7 1 Pb 



7 dp c \ _j , sin(/x' & 6) 
= 7 cos{/j c c)pb- 



(3.52) 



(3.53) 



2pl /2 dr 

Starting with / ~ and p c c ~ ^K c in a given cycle of classical phase, in the consecutive classical 
cycle across the big bounce, rather than constant, / turns out to be a widely time- varying function 
given as 



'PbPc +6 



7 2 A 



with some constant 5 to be determined. This follows 

COs(jj,' c c) ~ COS ±7T 



V ( R c + 8) 



PbPc 

cos{p! h b) w cos ( ±tt + (K' b + 5) 



(3.54) 

(3.55) 
(3.56) 



15 By contrast, in the //-scheme phenomenological dynamics of the Bianchi I model, / are indeed always constant (/ pa or / 2K./3) in 
classical regimes |10| . 
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and 



-i MKc) 

1 fc -, 



' PbPc 

'p 2 b Pc 

7 2 A 



sin ±7r + 



sin ±7r + 



PbPc 



PbPc 



-{K c + 8), 
iK + 5) 



(3.57) 
(3.58) 



provided the classical limit is held: p\p c ^> A r y 2 (K c + 5) 2 ~ Aj 2 (K b t ± + 5) 2 . Consequently, by (3.52) 
and (3.53), we obtain that in the consecutive cycle of classical evolution 



new effective K c 
new effective K b 



K c or K c 
K h or K h 



K c + 5, 
K b + 5, 



(3.59) 
(3.60) 



where K c (resp. K c ) and K b (resp. K b ) denote the new effective K c and K b in the next (resp. 
previous) classical cycle across the big bounce. Note that (3.54) indeed gives a new constant effective 



As in the /z-scheme analysis (see footnote 12), we assume that the bounce in p b take s plac e in 
the -fTfci-asymptotic phases (i.e. p\ <C K 2 + KT). Furthermore, notice that (2.31) and (3.49) are 



formally identical. As a result, once the universe enters the iT^-asymptotic phase, K b remains 
almost constant (K' b « K b ±) even when quantum corrections take effect later. (In the bouncing 
period, quantum effect varies / dramatically but modifies K' b only slightly.) 16 Exploiting this fact 



K 2 + Kl + 5, 



and using (3.59) and (3.60), we have 

K b _(K c ) = K b+ (K c )+5 

& -{K c + 5) + ^(K c + 5) 2 + K 2 = -K c - 
K b+ (K c ) = K b _(K c ) + 5 

& -(K c + 5)-^(K c + 5) 2 + K 2 = -K c +^K 2 + K 2 +5, 

which yields 5 = 2K c /3 ± A(K 2 C + K 2 ) 1 ' 2 /?, and gives the new effective K c : 



K r 



K r 



5 -K c + ^jK 2 + K 2 



(3.61) 
(3.62) 

(3.63) 
(3.64) 



Note that IQ = K c = K c . 

In summary, the effective K c in one cycle of classical regime is shifted to K c in the next classical 
cycle across the putative big crunch and to K c in the previous classical cycle across the putative big 



bang. Schematically, we have [cf. (3.32)] 



K r 



h. big crunch. big crunch - — ' big i 

■> K c , K c . K c — 



■* K r 



(3.65) 



16 Do not confuse K! with the effective K^. The former remains constant through the big bounce while the latter is offset by 8 as I 3.60 i 



suggests. 
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and [cf. gjjl] 



K 2 r + K 



big crunch 
2 ► 



K±.(K e ) = -K e + JK* + Kl, 



K b+ (K C ) = -K c -JK* + Kl ; 



big crunch 



K b -{K C 



(3.66) 



To find out the exact condition for the occurrence of big bounces, by substituting (3.47) and (3.48) 



into (3.41) and (3.43) and regarding cj) as the emergent time via (3.39), we study the differential 



equations: 



^dpc 2V4^G p 2 bPc 
p c d<t> V7 2 A 



cos 



' 7 2 A 

vlvc 



{K c + f) sin 



vlvc 



(K + f) 



(3.67) 



1 dp b y/AirG \ p\p c 

p~ b l4 ; KTy^fA 



cos 



^ (*a + f) 

PbPc 



X 



sin 



t£(K' b + f)) +s in(Jt£(K c + f) 



PbPc 



P b Pc 



(3.68) 



These are the modifications of the classical counterparts (2.37) and (2.38). [Also note that, as 
in the classical and /i-scheme dynamics, the flipping K c — > —K c t ogeth er wi th K' b ((p), f(<ft) — > 
—K' b (—(f)), —f(—<p) gives rise to the time reversal according to (3.67) and (3.68).] 



Similar to the case of (3.20) in the /z-scheme, p c gets bounced once the "cos(- • • )" term in (3.67) 



flips signs. This happens when 



cos 



IfA 

PbPc 



(K c + f) =0 



K c + f 



K PbPc 



2 V 7 2 A' 



(3.69) 



Assuming p b also gets bounced roughly around the same moment, 17 at which (3.69) is satisfied, we 
have the approximation: 



sml,/ lf (/?;+/) ) =si,i !2 



y± (K> - K c ) 

P b Pc 



COS 



PbPc 



ZPbPc 



(3.70) 



Taking (3.69) and (3.70) into (3.50), we have 

7 2 A 



} PbPc 



i - 



2p 2 bPl 



(K' b -K c 



+ 



PbPc 

7 2 A 



1 - 



7 2 A 

ZPbPc 



-i 2 



{K'b ~ K c 



+pi + 



(3-71) 



This is because I 3.67 i and (3.68 I are coupled through pw,. We will see that this is indeed the case in the numerical solution. 
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which, provided that ^> pb,p c ^ 7 2 A when the bounce occurs, 18 leads to the condition for the 
occurrence of the bounce: 

PbPc _ 1 

7 2 A ~ 6 

=: F(K c ,K' b ). (3.72) 

Since the Taylor series of cos a; = 1 — x 2 /2 + ■ ■ ■ converges very rapidly, the approximation made 
above is fairly accurate if 



2{K' b - K c f + K}+ J(2(K' b - K c y + K})' - 3(K' b - K c 



x 



' 7 2 A 

vlvc 



\K',-K r \ < 7T 



P 2 bPc 



-y 2 A - 



7T- 



(3.73) 



at the epoch when the bounce takes place. 19 

Similarly, p b gets bounced once the "cos(- ■ ■ 
argument above, we conclude that the big bounce of Pb happens when 

7 2 A , r>/ 



term in (3.68) flips signs. Following the same 



Ki 



which leads to 



,PbPc 

' 7 2 A 



PbPc 

7 2 A 



2p 2 b p, 



+ 



PbPc 

7 2 A 



+ Pt + --- 



\{K' b -K c f 



3 



G(K c ,K' b ) 



(3.74) 



(3.75) 



As mentioned earli er, we have assum ed that the big bounces take place in the -fT&-|--asymptotic 
phases; therefore, in (3.72) and (3.75), we can replace K' b with K b ±. To summarize, we conclude 

approaches the critical values 



that the big bounces of p c and pb take place when the matter density f 
p c ±,crit and p&±,crit, respectively, given by the Planckian density ppi := 
factors: 



Pc±, crit 



F(K c ,K b 
Kl 



Ppi, 



(87rC?7 A) times numerical 



(3.76) 



Pfcfc.c 



Ppi, 



G(K c ,K b± )^ (3 ' 77) 

where "+" is for the bounce which resolves the big crunch singularity and "— " is for the bounce 
which resolves the big bang singularity. 
The differential equations (3.39)-(3.43) can be solved numerically for given initial conditions. 20 



The numerical solution is shown in Fig. [3 The occurrence of big bounces is indicated by the matter 
density p^. The bounces of p b and p c take place when is close to p6±, C rit and p c ±,crit- Contrary 
to the /2-scheme, the epochs of bounces in p b and p c are roughly around the same time. Also notice 
that both p b and p c oscillate more and more rapidly toward the future and past; eventually pb 
grows very huge (and p^ subsides) while p c descends into the deep Planck regime, in which the 
quantum fluctuations becomes significant (in particular the cotriad component uj c = L^Jg xx grows 
huge and the quantum corrections on it have to be taken into account) and therefore the analysis of 
phenomenological dynamics can no longer be trusted. 21 This can be understood by the fact that the 



We will see that this is true until p c eventually descends into deep Planck regime in the far future and pase. (Cf. also see footnote [l2] 
for the case of the £Z-scheme.) 

Again, this can be justified in the end (cf. footnote fl2l . 

Unlike the case of the /I-scheme (see footnotc |l4| , the numerical method for the /2'-scheme phenomenological dynamics encounters no 
problems until the solution of p c descends into the deep Planck regime. 

As pi, and p c oscillate too fast, the numerical method also fails to give an accurate solution. 
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absolute value of the effective K c in the classical cycles becomes larger and larger toward the future 
and past and as a result the semiclassicality is less and less established. 

Contrary to the p-scheme dynamics, in which the directional densities g b and g c are the indications 
of bounces, it is the ordinary matter density p^ that signals the occurrence of bounces in the p'- 
scheme. Unlike g b , g c , the quantity of ps is independent of the choice of X since 



V 



2a2 



327r 2 p 2 p c V 2 



(3.78) 



Furthermore, (3.44) and (3.45) imply that the quantities p! b b,pf c c dep end o nl y on Pb/Pb, Pc/Pc and 
thus are independent of X (recall pb oc L, p c oc L°). Consequently, (3.47), (3.48) and ( 3.50) tel l 
us: K a , K' b , f and all scale as oc L. Therefore, the phenomenological dynamics given by (3.67) 



and (3.68) is completely independent of the choice of X as is the classical dynamics. In particular, 
the choice of X has no effect on the numerical factors K^/ F(K C , K b ±), K^/G(K C , K b ±) appearing 
in (3.76) and (3.77). This is a desirable feature that the p-scheme does not have. [However, if we 
further impose the quantum corrections on the eigenvalue of the cotriad operator u> c , this invariance 
is broken again.] 

Even though the p'-scheme is independent of X, in case when the finite sized shell X x S 2 has a 
global meaning, the bounce conditions (p^ ~ p c ±,crit, P&±,crit) can be understood as: The volume of 
X x S 2 (i.e. V = ^npb\fp~cj gets bounced when it undergoes the Planck regime (times numerical 
values F X I 2 {K C ^ Kb±)/K ( p, G 1//2 (i^ c , K b ±)/K^) measured by the reference of the momentum p^. 



IV. SCALING SYMMETRY AND RELATIONAL MEASUREMENTS 

With LQC discreteness corrections, the phenomenological dynamics (both in the p- and p'- 
schemes) shows that both the big bang and big crunch singularities are resolved and replaced by 
the big bounces. The occurrence of bounces is indicated by the directional densities gb, g c in the 
p-scheme whereas it is signaled by the matter density p^ in the p'-scheme; the bounces take place 
when gb, g c in the p-scheme or p^ in the p'-scheme approaches the critical values. 

It has also been noted that the classical dynamics and the phenomenological dynamics in the 
p'-scheme are both completely independent of the choice of the finite sized interval X, while the 
phenomenological dynamics in the p-scheme reacts to the physical size of X. This can be rephrased 
in terms of the scaling symmetry; 22 that is, the classical dynamics and p'-scheme phenomenological 
dynamics are invariant under the following scaling: 

Pb, Pc — > Ipb, Pc, 

b, c — > b, Ic, 

p^, Kf — > Ipcj,, IK^, 

K c — > IK C . (4.1) 

[Note that the scaling for K c should be accompanied by the same scaling on K b in classical dynamics 
and on K' b as well / in the p'-scheme; that is K b , K' b , f — ► lK b , lK' b , If.] On the other hand, the 
p-scheme does not respect this scaling. 23 



22 A dynamical system is said to be invariant under a certain scaling if for a given solution (pi,(r), p c {t), b(r), c(t), <f>(r) and p^) to 
the dynamics, the rescaled functions also satisfy the equations of motion (i.e. Hamilton's equations and vanishing of Hamiltonian 



constraint). For t he cla ssical dyn amics, the equations to be satisfied are ! |2.21| H |2~27l l; for the ^t-scheme, pli|)-| |3.1l| ); and for the 
/t'-scheme, | |3.38| -| (3^43| and ( |3.50 I ■ 
23 Also note that the symmetry involving <j> * ^<j> as stated in Equation (4.2) of I10| is a mistake and should be dismissed. 
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FIG. 3: Solution in the //-scheme phenomenological dynamics. Same initial condition as given in Fig. [j] (and 7=1). 
(a) pb(4>) and Pc(</>) = gnn(<p)- The bounces of pb and p c occur around the same moments, pb and p c both oscillate (more and 
more rapidly) toward the future and past: pb grows up while p c descends into deep Planck regime, at which the semiclassical 
description can no longer be trusted, (b) g xx {4 1 )- ( c ) P4> (</>)> which signals the occurrence of big bounces. The critical values 
given by (3.761 and (3.771 with different effective K c are indicated by dashed lines. The epochs of big bounces are very close 



to the moments when the curve intersects one of the dashed lines, (d) cos(fi' b b) and cos(//c), fairly close to each other, 
(e) f(4>). f ~ in the classical cycle around (j> = (j>o but / yields the time- varying signature given by (3.54| in the adjacent 



classical cycles, (f) K' b ((j)), which remains almost constant across big bounces, (g) Effective K c with the constants given by 
(3.651 indicated by dashed lines, (h) Effective K b with the asymptotic values given by (3.661 indicated by dashed lines. 
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The fact that the symmetry in (4.1 ) only scales pb and c but leaves p c and b fixed seems to spoil the 
idea that length/area/ volume is measurable only if the line/surface/bulk is coupled with the material 
reference as suggested in [2H 122] and affirmed in [TUj for the Bianchi I model. This is because in the 
Kantowski-Sachs spacetime the area = Tignn (contrary to A^) has no ambiguity even in the 
absence of matter content; in a sense, A^ is measurable with reference to the spherical curvature of 
S 2 and thus the reference to p$ is unnecessary. In the l£&±-asymptotic phases, where the spherical 
curvature is negligible (compared to the anisotropy), we expect the same scaling symmetry as of 
the Bianchi I model (see Equation (4.1) in [10] and recall Remark [l] in Sec. II B). That is, if we 



ignore ^p 2 in the bracket of (2.20) and the corresponding term in the effective Hamiltonian, the 
classical dynamics and /i'-scheme dynamics in the i£y_-asymptotic phases will approximately respect 
the symmetry of scaling under: 



Pb, Pc — ► X b X c p b , 



b, c — > Xbb, X c c 
\ 2 



K<f> — > XlXcPtf,, X\X C K^ 



K c — > X 2 b X c K c . (4.2) 

The /2-scheme dynamics respects this approximate symmetry as well if we set Xb = X c . 
Additionally, the classical dynamics also admits the symmetries given by 

T ► 7]T, 

7 — ► £7, 
Pb, Pc — ► V 2 Pb, V 2 Pc 
b, c — ► £b, fc, 

P<t>, K <t> — * V 2 P<t>, tf K <i>, 

K c — > V 2 K C . (4.3) 

The scaling symmetry regarding 7 — > £7 is expected, since the Barbero-Immirzi parameter 7 
has no effect on the classical dynamics. The scaling symmetry regarding r — ► rjr is also easy 
to understand, since there is no temporal scale introduced in the Hamiltonian. 24 However, very 
surprisingly, the scaling symmetry involving r — ► i]T is violated for both the /i-scheme and /i'- 
scheme phenomenological dynamics. Curiously, this symmetry is restored if r — ► ryr is accompanied 
by 7 — ► £7 and one extra scaling is also imposed at the same time: 

A — ► rVA- (4.4) 

This intriguing observation seems to suggest, albeit speculatively, that in the context of quantum 
gravity the fundamental scale (area gap) in spatial geometry gives rise to a temporal scale via 
the nonlocality of quantum gravity (i.e., using holonomies) and the Barbero-Immirzi parameter 7 
somehow plays the role bridging the scalings in time and space. [This reminds us that, in LQG, the 
precise value of the area gap A is proportional to 7, and 7 is also the parameter which relates 
the intrinsic geometry (encoded by spin connection r*) with the extrinsic curvature (K a % ) via 

iC = P.-7tf.*-] 



24 For the Bianchi I cosmology studied in |10| , a different scaling pj > pj with cj > r)~ 1 cj is chosen to respect the symmetry regarding 

r > ijt. This alternative scaling does not work in the case of Kantowski-Sachs spacetime, since it violat es the Hamiltonian constraint 

( 2.20 k That is to say, the presence of the spatial curvature [i.e., the 7 2 p 2 term in the bracket in 1 2.20 1] ties the temporal scale with 

the spatial scale; as a result, only the scaling pt, p c ► J ) 2 P6i ^ 2 Pc (with gives the spatial direction the same scaling as in the temporal 

direction) with b, c > b, c preserves the symmetry. 
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Meanwhile, related to the above observations, the physical meaning of the directional factors qj 
and matter density p^ can be interpreted as the (inverse of) area and volume scales, measured by the 
reference of the matter content. In this regard, we may say that the big bounces take place when one 
of the areas A^, A X( p (in the /i-scheme) or the volume V (in the /i'-scheme) undergoes the Planck 
regime (up to a numerical factor) measured by the reference of the matter momentum p^. It is then 
tempting to regard not only as the "internal clock" (emergent time) but also p^ as the "internal 
rod" — namely, the measurement of both temporal and spatial geometries makes sense only in the 
presence of matter content. 

The above observation for the scaling symmetries draws a close parallel to those in the Bianchi 
I model [10J and seems again to support the ideas of the relational interpretation of quantum me- 
chanics with real rods and clocks such as studied in [23] (see also [2TJ [22]) with the caveat that the 
nonvanishing spherical curvature can make the reference to matter content unnecessary. The caveat 
however does not fail the relational interpretation immediately; rather, it suggests that we should 
put the spatial curvature on the equal footing as the matter content. 25 More precisely, apart from 
energy density of matters to be the metric reference, we should also take into account the energy 
densities of both curvature and anisotropic shear. 26 

Unfortunately, all the scaling symmetries break down in the detailed construction of LQC with 
the /i-scheme (the strategy to construct the fundamental theory of LQC in the /i'-scheme is still not 
clear) even for the isotropic model (where the p- and /i'-schemes are identical). The fundamental 
LQC only respects the scaling symmetries at the leading order. This is due to the fact that the 
quantum evolution in the fundamental LQC is governed by a difference equation, in which the step 
size of difference introduces an additional scale in the deep Planck regime (see [I] for the isotropic 
model and [7] for Bianchi I model). In fact, already in the level of phenomenological dynamics, the 
scaling symmetries are violated if we further take into account the LQC corrections on the cotriad 
component u c . For the fundamental theory of LQC, if we take the aforementioned symmetries 
seriously, we might be able to revise the detailed construction in the spirit of relational quantum 
theory such that the step size in the difference equation scales adaptively by the reference of the 
matter content. 

V. SUMMARY AND DISCUSSION 

To summarize, we list the important facts for the classical dynamics, /i-scheme and /i'-scheme 
phenomenological dynamics in Table [Tj In the following, the main results are restated and their 
implications are discussed. 

With the LQC discreteness corrections, the phenomenological dynamics shows that the classical 
singularities (both big bang and big crunch) are resolved and replaced by the big bounces. In the 
/i-scheme, it is the directional densities Qb and g c that signal the occurrence of big bounces when Qb 
and q c approach the critical values ^6±, crit an d £> C) crit respectively. In the /i'-scheme, the indication 
of big bounces is the matter density p^ and the big bounces take place around the moments when 
Ptf, is close to the critical values pb±, crit and p c ± >cr it- 

The detailed evolution in the /i-scheme shows that the equations of motion (in terms of emergent 
time (p) for pt, and Kb are decoupled from p c and c (the dependence on p c and c is only through 
the constant K c ). As a result, the bouncing scenario of pb is unaffected by the varying of p c and 



This reminds us that in the FRW model the curvature term for k = ±1 can be regarded as "matter" with state parameter w = —1/3. 
It has been shown in Appendix B of [9] that the anisotropic shear behaves as a kind of anisotropic matter; in particular, the directional 
densities can be considered as the "energy density carried from the classical anisotropic shear portioned to the specific direction" . 
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Classical dynamics 


Phenomenology in /i-scheme 


Phenomenology in ^'-scheme 


P4 , = V47TG- 1 A = W<f> 


P4, = V4-KG- 1 = Vcf> 


p = V47TG- 1 A = V<f> 


2%p b = 2TvL^/g xx gnn = A^ = A x e 
np c = 7r(/nn = A<^ 


2np b = 2-KLy/g xx gaa = A xl j, = A x e 
■npc = Tignn = Agj, 


2np b = 2-kL^J g xx gan = = A x s 
7rp c = Ttgna = A^ 




sin(/i b b) 1 7 dp c 


sin(pj,6) _ l 7 dp e 




_ pj/ 2 dr 2p 3/2 dr 


fib ~ cos(m c c) 2p l /2 dr 
sin(p, c c) _ 1 7 dpp 


A' b cos(££c) 2p y 2 dT 
sin(/i^.c) i 7 dp p 


Me ~ co${[i b b) p l/2 dr 

1 7Pb dpc 


M c _ cosfp^b) p l/2 dr 

1 7Pb dp e 
cos(p^c) 2p 3 ^ 2 df 


p c c = 7A C 


p c sin i^ c) = 7 A C 


p c c = 7 [Ac + /(</»)] 
p 6 6 = 7[A^(<A) + /(0)] 


Kl = 2K b K c + Kl+pl 


= 2K b K c + K 2 b + p 2 b 


P^Pc 

2sin (\/l? ( ^ +/) ) sin (\/l? ( ^ +/) ) 
+ sin2 (v / S^ +/ )) + i? 


K 4> dK b _ 2 




P^ = ^lf V^cos(/|A ( ^ + /) 
x sin (\/S^ + / )) 

x [ sin (\/s ( ^ +/) ) 
+ sin (v / S^ ( ^ +/) ) 


) 


p c , p b -» 

toward both big bang and 
big crunch singularities. 


p c bounces whenever 

£ c — 32,r 2 p3 - KjPPh 

p b bounces whenever 

4 ~ K l n 

Sb ~ 3 2^ p 3 « Jf2 ± />P1, 

P6 bounces periodically; 

p c bounces only a few times. 


p c bounces at the moment when 

4 ~ K l 

P<t> — 32^2 p 2 pc - F(Kc ,K b ±)PPU 

p b bounces at the moment when 

P<t> 32^p'^p c G(K c ,K b± )P pl > 

Pb, p c bounce roughly around 
the same moments. 


No big bounce; 
A c fixed 


Big bounces bridge classical cycles 
with K c and — A c and give rise 
to "meta-classical" phases. 


Big bounces bridge classical cycles 
with varying effective A c : 

• • • Ac t± Kc t± Ac t± Kc t± Kc ■■■ 


Symmetry of scaling: 
r > yyr 

7 — > C7 
Pb, p c — > /r? 2 p 6 , ?7 2 p c 
6, c — > £&, Z£c 
p<f>, K^> — > lri 2 p$, IrfKj, 
K c — > lifK c 


Symmetry of scaling: 
T > rjT 

7 — > 

Pb, Pc > tfpb, J? 2 Pc 

b, c — > £6, £c 
P0, A^ — > 7) 2 ^, j^A^, 

A c > ?) 2 A"c 

a — r Va 


Symmetry of scaling: 

t > rjT 

7 — > 

Pb, Pc — > Irfpb, V 2 Pc 
b, c — > £b, l£c 
p<p, A^ — > Irfp^, lr/ 2 Kj, 
(Kc + f)^l V 2 (Kc + f) 

A — ► tVa 


Approximate symmetry in 
A(,±-asymptotic phases: 
Pb, Pc — > Xb^cPb, X\p c 
b, c — > X b b, X c c 
P<P, K<t> — > X\\ c p$, X b X c K$ 
K c — ► A 2 A C A C 


Approximate symmetry in 
A(,±-asymptotic phases: 
p b , p c — > A 2 p b , A 2 p c 
6, c — > Xb, Ac 
p^, A — > A 3 p^,, A 3 A^ 
Ac — > A 3 Ac 


Approximate symmetry in 
At±-asymptotic phases: 
Pb, Pc — > AbAcPf,, A 2 pc 
b, c — > X b b, X c c 
P4>, Kj, — > A 2 A c p0, A 2 A C A^ 
(Ac + /)^A 2 A c (Ac + /) 



TABLE I: Summary of the classical dynamics, ^-scheme and //-scheme phenomenological dynamics. 
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perfectly periodic. On the other hand, p c only bounces a few times and grows up toward infinity in 
the far future and past. By contrast, in the /i'-scheme, equations of motion for pb and p c are closely 
coupled through p^ and thus pb and p c bounce roughly around the same moment. 

The big bounces bridge different "cycles" of classical solutions. In the /2-scheme, the classical 
solutions with (effective) constants K c and — K c are bridged by the big bounces. For some periods 
of time, the phenomenological dynamics also yields "meta-classical" phases, which are absent in the 
ordinary classical evolution. (Further investigation is needed to know whether the occurrence of 
meta-classical phases indicates wrong semiclassical behavior.) On the other hand, in the /i'-scheme, 
starting with the constant K c in a given classical cycle, the evolution ends up with the new effective 
constant K c or K c in the adjacent classical cycle across the big bounce. Furthermore, the bouncing 
behaviors of pb and p c oscillate more and more rapidly toward the future and past; eventually pb grows 
very huge (and p^ subsides) while p c descends into the deep Planck regime, where the semiclassical 
analysis of phenomenological dynamics can no longer be trusted and the quantum corrections on 
the cotriad component u c become important. 

In regard to the finite sized interval X chosen to make sense of the Hamiltonian formalism, the 
phenomenological dynamics in the /2-scheme depends on the choice of X, and hence reacts to the 
macroscopic scales introduced by the boundary condition. (In terms of symmetry, it is said that 
the /Z-scheme has no such scaling symmetry respected by the /i'-scheme and classical dynamics.) 
The phenomenological dynamics in the /i'-scheme, in contrast, is completely independent of X as 
is the classical dynamics. In case that the physical size of X has a global meaning (such as in the 
compactified Kantowski-Sachs model or in the lattice refining model of [20J ) , the condition for the 
bounce occurrence can be rephrased: In the ////i'-scheme (respectively), the physical area/volume of 
the surfaces/volume of A x q, or V gets bounced when it undergoes the Planck regime (times a 
numerical value) measured by the reference of the momentum p^. 

While the /i'-scheme has the advantage that its phenomenological dynamics is independent of X, 
the fundamental theory of LQC based on the /i'-scheme is difficult to construct. Both the p- and 
/Z'-schemes have desirable merits and it is still disputable which one (or yet another possibility) is 
more faithful to implement the underlying physics of quantum geometry. 

In addition to the symmetry related to the choice of X, both schemes admit additional symmetries 
of scaling, which are reminiscent of the relational interpretation of quantum mechanics, featuring the 
ideas of real rods and clocks. Furthermore, the symmetry involving the Barbero-Immirzi parameter is 
suggestive that the fundamental scale (area gap) in spatial geometry may give rise to a fundamental 
scale in temporal measurement. These symmetries however break down in the construction for the 
fundamental theory of LQC. 

Most results obtained in [TU] for the Bianchi I model are analogously affirmed at the level of 
phenomenological dynamics for the Kantowski-Sachs spacetime and the close parallels between these 
two cases are well established. (Also see Remark [I] inSec. II B ) Additionally, new features also arise 
due to the presence of the spherical curvature: The bouncing scenario exhibits (semi)-cyclic patterns; 
and in the /i'-scheme, p c eventually descends into the deep Planck regime, whereby the validity of the 
phenomenological dynamics could be questioned and more sophisticated treatment may be required 
to faithfully convey the quantum geometry of the full theory of LQG. At this stage, it is not clear 
what happens exactly when p c reaches Planck regime. 

Meanwhile, it is noteworthy that the Kantowski-Sachs spacetime describes the interior of the 
Schwarzschild black hole (see Remark [2] in Sec. II B). Although, by introducing the scalar field, 
we dismiss the issues for black holes and instead study a self-contained cosmological model, it is 
still instructive to compare the phenomenological dynamics of the Schwarzschild interior studied 
in [15] with the results obtained in this paper. It has been shown in [15] that in both the p Q - 
and /2-schemes (referred to as "constant 5 Hamiltonian" and "alternative quantum Hamiltonian" 
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in [H]), the phenomenological dynamics bridges a classical black hole with a white hole through a 
bounce, whereas in the /i'-scheme phenomenological dynamics (referred to as "improved quantum 
Hamiltonian" in [IS]), Pc oscillates and eventually lands on a constant in deep Planck regime. 27 This 
dichotomy is analogous to the qualitative difference in the bouncing cosmological scenarios between 
the pr and /l'-schemes observed in this paper. This similarity also suggests that, by exploiting the 
strategy we used in this paper (in particular, to identify the constant K c and the functions Kb, K, 
K'), we could be able to reproduce the results of [15] for the black hole interior in greater detail (such 
as to pinpoint the bounce occurrence condition and the varying of the effective K c ); the quantum 
corrections on the inner side of horizon could also be studied. 

Furthermore, as this paper focuses specifically on the model with a massless scalar field, it should 
be straightforward (with necessary approximation) to extend the results to the models with inclusion 
of generic matters. As studied in [5] for the Bianchi I model, it is anticipated that there would be 
a competition among the matter density, anisotropy and spherical curvature to be the indication of 
the occurrence of big bounces. Studying the loop quantum geometry of Kantowski-Sachs spacetime 
with generic matters would further support or oppose our observations and help us to understand 
them in a broader context. 
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APPENDIX A: PHENOMENOLOGICAL DYNAMICS IN THE ^-SCHEME 

One of the virtues of the improved strategy (ft- or /i'-scheme) in both the isotropic and Bianchi I 
models is to fix the serious drawback in the old precursor strategy (/i Q -scheme) that the critical value 
of directional densities qi (in the /i-scheme) or of matter density (in the //'-scheme) at which the 
bounce occurs can be made arbitrarily small by increasing the momentum p^, thereby leading to 
wrong semiclassical behavior. 

Having learned from the isotropic and Bianchi I cases, we expect that the critical values of g c , q 
and at which the bounces occur can be made arbitrarily small by increasing the momentum p^ in 
the /io-scheme but are independent of p^ in the /2- or /i'-scheme. The latter is what has been shown 
in the main text of this paper. 28 For comparison, the phenomenological dynamics in the /i -scheme 
is presented here. 

In the phenomenological theory of the /i -scheme, we take the prescription to replace c and b with 
sin(/i°c)//x° and sin(/i£6)//z£ by introducing the fixed numbers /i° and p° h for discreteness. Analogous 



Therefore, it is claimed in [15j that the /i-scheme phenomenological dynamics extends a classical Schwarzschild black hole to a patch of 
a nonsingular charged Nariai universe, which gives constant p c . However, a closer look suggests that the extended part is not a patch of 
the classical Nariai universe but instead represents the quantum universe which formally exhibits Nariai type metric, as the asymptotic 
constant value for p c is in the deep Planck regime (< ^j). For this reason, some of the claims in I15| for the /t'-scheme may require 
further investigation. 

In the /i-scheme, the critical values £> c ,crit and £>j,±. cr it depend on p^ only through the ratios K^/K^ and K^/K^; In the /i'-scheme, 
Pc±,crit and p b ± iCr it depends on p^ only through K?/F(K C , K b ±) and K^/G(K C , K b± ). 
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to (3.2), we have the effective (rescaled) Hamiltonian constraint: 

f 



2G 7 



sin(>£6) sin(^c) 
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Again, the equations of motion are given by the Hamiltonian constraint H'^ o = and Hamilton's 



equations: 
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(A2) 
(A3) 
(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 



These are exactly the same as (3.14)-(3.17) except that the discreteness variables p, c , fib are now 
replaced by \x° c and ^tjj. 

Therefore, exploiting the close resemblance between the /2-scheme and /i D -scheme, we can readily 
repeat the calculation we did in Sec. Ill A and obtain the differential equations [cf. ( |3.19 )-( p.21[ )]: 
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which give the bouncing solutions similar to those given in the /i-scheme phenomenological dynamics 
except that the critical values at which the big bounce takes place are given by 
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which can be made arbitrarily small by increasing the value of p^. As a result, the p -scheme gives 
wrong semiclassical behavior and should be improved by the p> or p'-scheme to fix the problem. 



APPENDIX B: EFFECTIVE HAMILTONIAN IN THE /i-SCHEMES 



In this appendix, starting from the Hamiltonian constraint of LQG, we derive the gravitational 
part of the Hamiltonian with Kantowski-Sachs symmetry and give a heuristic argument for the 
prescription given by (3.1). The motivations for both the p- and p'-schemes are addressed in detail. 
The advantages and drawbacks of both schemes are also remarked. 

The gravitational part of the classical Hamiltonian constraint in the full (unreduced) theory is 
given by 



H, 



grav 



^ J d'xNe- 1 {e^ k Fi b E a 3 E b k - 2(1 + 7 2 ) W E\e\} , (Bl) 



where e := | det E] 1 / 2 = ■ s /q. 

With the Kantowski-Sachs symmetry, the connection potential is given by (2.6), which leads to 
the field strength: 



F = -F ab dx a A dx = dA + A A A 

= (b 2 - l)r 3 sin 9 d9 A dcj) + 2bcr 2 sin 6 d(j) A dx - 2bcr 1 dx A d9. 



(B2) 



On the other hand, the densitized triad given by (2.7) gives e = y/q = pb\^p\sm9 and the corre- 
sponding cotriad: 



UJ = UJ % „T;dx a 



Pb 



=r%dx + \/PcT~2d9 — \/p c T\ sin 



'Vc 



via uo l a E j = y/qdj. The compatibility relation duo + V A uo = then yields the spin connection: 

r = T3 cos 9d<p. 
Consequently, the extrinsic curvature K is given by 

•yK := A — T = cr 3 dx + br 2 d9 — br\ sin 9 dcj), 

which follows 

y 2 K A K = b 2 T3 sin 9 d9 A dcj) + b&T2 sin 9 dcj) A dx — bc.T\dx A d6. 
Therefore, we have 
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€i jk F^ b E a jE b k = (b 2 - l)p 2 b sin 2 9 + 2bdp b p c sin 2 9 
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and 



1 2 K { :K h ^E a , i E\ = H 2 pl sin 2 9 + bcp b p c sin 2 9. 



(B8) 



Putting (B7) and (B8) into (Bl) and restricting the integral to the finite sized shell Z x S as 
prescribed in (2.8), we then have the gravitational part of the classical Hamiltonian in terms of the 

N 



reduced variables: 
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which is the H gISLV given in (2.16). 

When the quantization is performed in the context of LQC, there are two loop quantum corrections. 
The first is the modification on the cotriad operator u> c , which is negligible and ignored in this paper. 
The second is due to the fact that the connections A a l (or b, c) do not exist and should be replaced 
by holonomies (or exponentials of b, c). 

Following the standard techniques in gauge theories, the curvature component F % ah can be expressed 
in terms of holonomies (i.e. Wilson loops). Given a small surface a center in x, Stokes' theorem 
allows us to write the curvature component as 



r i F t cd dx c Adx" 



~ab 



Vexp 
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where da is the boundary loop of a and e^ b = f dx a A dx b is the coordinate area of a projected in 
"a6-direction" . This is a good approximation provided e" b is small enough and in fact it becomes 
exact in the continuous limit e" b — > 0. 

With the Kantowski-Sachs symmetry, we choose a to be small rectangular surfaces D^, and 
□ a normal to the vectors d x , dg and respectively. The coordinate lengths of the edges of a along 
the directions d x , dg and are denoted as p, x L, fig and jx^/ sin# with the discreteness parameters 
introduced. (See Fig. ma).) We then read off from (BIO) that 
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with h x ^ x \ and h^' being the holonomies along the individual edge of a: 
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Essentially, this is to replace the components of the connection A in the form of (2.6) with the 
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where the extra factor 2 in h x , h e , and /i^ is adopted to be consistent with (Bll) 
According to (Bl), (B2) and (B6), the relevant components of F and K A K appearing in -n gr av are 
6 2 r 3 sine, bcr 2 sine and —bcrx, which now can be expressed in terms of holonomies as 
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Equivalently, to take into account the discreteness corrections of LQG, we make the following pre- 
scription: 
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The continuous limit with p, x ,flo, p>cj> recovers the classical Hamiltonian constraint in (B9). 
However, the very feature of LGC is that the continuous limit does not exist and the failure of the 
limit to exist is intimately related with the underlying quantum geometry of LQG, where eigenvalues 
of the area operator are discrete and has an area gap A. In LQC, to implement the discreteness as 
imprint from the full theory of LQG, we have to set the discreteness parameters p, x , fig and to 
be finite. There are many possibilities to fix the discreteness parameters, but as in the Bianchi I 
model (see Appendix B of [10J), two well-motivated strategies (/2- and /i'-schemes) are of particular 
interest and presented in the following. 

The first strategy to impose the discreteness of LQG is to set the physical areas of D^, and 
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/sing 




FIG. 4: (a) The surface in pink is □0^, the physical area of which is to be shrunk to A in the fj! -scheme, (b) The surfaces in 
blue are Elg and H^, the physical areas of which are to be shrunk to A in the p-scheme. 



n x g to be A (depicted in Fig. [5](a)). That is 
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which, along with (B23), gives the "/i'-scheme" in (3.4). 

Instead of shrinking the areas of D^, and D x9 to A, the second strategy is to associate each 
edge of D^, D^. and H x q with an area and then shrink the associated areas to A. (See Fig. |4^b).) 
For instance, The edge of EJ^ in ^-direction is of coordinate length fig, with which, most naturally, 
we associate a rectangle E\g normal to ^-direction of coordinate lengths figL and fig/ sin8 on its 
edges. We then set the physical area of E\g (and similarly of and E\ x as well) to A. That is 
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which, along with (B23), gives the "/i-scheme" in (3.3). 

[Note that in both pictures of the fi'- and /i-schemes, if we rather set the coordinate areas of 
□e^, • • • or E\g, ■ ■ ■ to be A, we end up with constant fib and fi c . This is the "/i -scheme" used in the 
precursor strategy. The wrong semiclassical behavior of the /i D -scheme partly originates from the 
problem that the coordinate area (both of D^, ■ ■ ■ and of E\g, ■ ■ • ) is not totally physically relevant.] 

In the full theory of LQG, when the Hamiltonian acts on an spin network state, it adds a new 
link of spin-1/2 and the coloring of the links on which the new link is attached is increased or 
decreased by 1/2. The resulting spin network state is equivalent to the original state superimposed 
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with a triangular loop of spin- 1/2 links. This triangular loop is essentially the Wilson loop discussed 
above. In the context of spin network states, the coloring of a link corresponds to the area of the 
surface which the link penetrates and the smallest coloring spin-1/2 gives rise to the area gap A. 
Therefore, it is in this sense that the /i-scheme is a more direct implementation of the underlying 
discreteness of quantum geometry than the /i'-scheme. (See FIG. 6 and the pertinent text in [TU] for 
more comments.) 

Furthermore, in the fundamental quantum theory of LQC, the /i-scheme has the important virtue 
that we can define the affine variables v b , v c via 
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such that Hamiltonian constraint of the fundamental quantum theory gives the evolution as a dif- 
ference equation in terms of Vb, v c and therefore the methodology used for the isotropic model [I] 
and Bianchi I model [7j should be easily applied. This strategy fails in the /i'-scheme since 
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and hence the corresponding affine variables do not exist. This makes it difficult to construct the 
fundamental quantum theory in the /i'-scheme. 



On the other hand, as studied in Sec. Ill, the /i'-scheme has the advantage over the /i-scheme that 
the phenomenological dynamics in the /i'-scheme is independent of the choice of X. The difference 
for this point between the fl- and /i'-schemes can be understood, heuristically but instructively, by 



estimating the quantities and fx c c with the classical formulae; that is, substituting (2.28) and 
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Since the quantities fibb and /i c c indicate how significant the quantum corrections are (quantum 
corrections are negligible if £Lbb,p> c c < 1 ), (B30) and (B31) tell that in the /i-scheme, the place at 



which the quantum effects become appreciable is tied up with not only t he "H ubble rates " for y/gnn 
and y/g xx but also the physical geometry of X x S 2 . On the other hand, (B32) and (B33) show that 
p! b b and Jx' c c are proportional only to the "Hubble rates" in the classical regime and t 
choice of X is irrelevant in the /i'-scheme. 



re different 



Finally, it should be noted that, in the language of the lattice refining model in [20] . the /i-scheme 
corresponds to the refinement pattern that the number of lattice vertices is proportional to the 
transverse area, reminiscent of the idea depicted in Fig. |4j^b) with the transverse surfaces shrunk to 
A. A stability analysis however suggests that the /i-scheme leads to an unstable difference equation 
of evolution in the fundamental theory of LQC and thus, in this regard, could be problematic as a 
good quantization scheme. 
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Both the p,- and /I'-schemes have desirable and undesirable features of their own. In order to 
understood them more deeply, in the main text we study both schemes and their ramifications at 
the level of phenomenological dynamics. 
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